\(^1\)Universidad Autónoma de Madrid. Departamento de Prehistoria y Arqueología, Campus de Cantoblanco, 28049 Madrid, Spain.
\(^2\)Institut Català de Paleoecologia Humana i Evolució Social (IPHES), Zona Educacional 4, Campus Sescelades URV (Edifici W3), 43007 Tarragona, Spain.
\(^3\)Àrea de Prehistoria, Universitat Rovira i Virgili (URV), Avinguda de Catalunya 35, 43002 Tarragona, Spain.
Abstract
Lithic artifacts are some of the most common and numerous remains from paleolithic archaeological sites. However, their deposition into the archaeological record can be the result of multiple episodes of discard which are
Because the numerous amount of lithic remains, a quick, flexible, low cost and effective method for identifying degrees of alteration in the surface of lithic implements is highly desirable. This would favor the analysis of complete lithic assemblages, inferring their overall integrity, identify different episodes or clusters of alteration, and even identify post-depositional alterations of different nature (such as rounding, dragging or wind erosion).
An experimental collection of flints and quartzites were subjected to sequential episodes of rounding. After each episode, photographs with a Dino-Lite Edge 3.0 AM73915MZT USB microscope were taken, allowing to obtain quantitative values of surface texture. Following this, the quantitative data of surface texture was employed as variables in Machine Learning models to determine degree of rounding and the most important variables for discrimination.
Results show, that despite being a continuum,
Key words: experimental archaeology; postdepositional alterations
Lithic artifacts constitute one of the most common remains from paleolithic archaeological sites. When analyzing a lithic assemblages, a key factor to consider is determining the degree of post-depositional alterations undergone by the assemblage since this results in behavioral and eco-cultural inferences. Methods for determining assemblage integrity usually focus on the assemblage as a hole. These methods often include spatial analysis of artifacts and analysis of fabrics to determine if water flow has resulted in a reorganization of the spatial distribution (Lenoble and Bertran, 2004; McPherron, 2005; McPherron, 2018; Petraglia and Potts, 1994; Schick, 1987). Also lithic size distribution analysis are realized to determine if postdepositonal processes have resulted in sorting of the materials (Bertran et al., 2012; Maíllo Fernández, 1998; Schick, 1987). Another option is to focus directly on individual lithic artifacts to determine the degree of alteration that they have undergone (Chu et al., 2015; Levi Sala, 1986). Although being more time consuming, this type of analysis provides a higher resolution, allowing to establish a graduation in the degree of postdepositional alterations, to combine its data with spatial analysis, dissecting archaeological episodes, or to detect several episodes of recycling.
Individually determining the degree of alteration undergone by a lithic artifact uses microscopic analysis to measure ridge width, visually analyze the surface to determine the existence and intensity of abrasion, and examine the edges to determine the presence of detachments coming from particle impact or dulling of the edges in more extreme cases (Bustos-Pérez et al., 2019; Chambers, 2003; Shackley, 1974). Of these measures, ridge width is the most clearly quantifiable which can result in an greater weight on the analyst’s decision. However, analysis from sequential experimentations shows that the longer the exposure time to sedimentary abrasion, greater the increase in the standard deviation of ridge width measurements (Bustos-Pérez et al., 2019). This introduces a degree of uncertainty when determining the degree of rounding undergone by a stone tool, since the ratio of increasing ridge width might be the result of general flake topography, ridge position among the flake, ridge angle, or general profile of the ridge.
Adding to the variability in flake width (which comes from flake and ridge morphology), it is also necessary to consider the effects of fluvial sedimentary processes which commonly affect lithic artifacts. Fluvial sedimentary processes are characterized by particle transport and these processes are not homogeneous. Lithic artifacts can form part of these process as another particle or remain static and be affected by the particles forming part of them. Commonly, three modes of particle transport are described in fluvial sedimentary contexts: rolling, sliding and saltation (Alhusban and Valyrakis, 2021), although the last one is not a common form of alteration in the case of stone tools (Schick, 1986). Additional to alterations coming from fluvial contexts, stone artifacts might also be affected by aeolian particle transport, usually resulting in wind erosion Stapert (1976).
Stone tools have an additional factor of variability added to the ones derived from particle transport and sedimentary contexts (fluvial or eolian). This additional factor comes from the degree of exposure and speed on which a stone tool enters the archaeological record (Petraglia and Potts, 1994; Schick, 1986; Schiffer, 1972). For example, it is expected that partially buried artifacts with water and sediment flowing above them will present modifications in the exposed surface, while edges and the burred surface will remain semi-intact (Petraglia and Potts, 1994). A different scenario happens when artifacts are transported by rolling in coarse sediments, resulting in abrasion of all surfaces, dulling of the edges, but at the same time, the impact from coarse particles might result in freshly detached surfaces and edges which undergo new modifications until the artifact enters the archaeological record (Harding et al., 1987; Petraglia and Potts, 1994).
Thus, obtaining quantitative data of surface modifications can complement values of ridge width, and help identify and interpret type of particle transport undergone by, or affecting a stone artifact. Ideally, quantitative values of surface abrasion would be extracted from each stone artifact, allowing to maximize information and interpretation of the formation process of an archaeological site. However, as previously mentioned, stone artifacts are some of the most common remains from Paleolithic archaeological sites. Thus, a versatile, fast an non-time consuming method is highly desirable since it would enable to analyze large quantities of lithic materials.
This research presents results from a sequential experimentation in which a set of experimentally knapped flakes have undergone rounding process in a tumbler. After each episode of rounding, the same areas are photographed, and quantitative measures of flake texture are employed to characterize the surface. Following the extraction of quantitative values, Machine Learning models are employed to determine the degree of separation between rounding times. Flake ridge width is also considered as an indicative measure of alteration undergone by flakes.
Three blocks of flint representing three different types (F1, F2 and F3) and three blocks of quartzite also representing three different types (Q1, Q2 and Q3) were experimentally knapped. Flint type F1 is characterize by a heterogeneous grey color surface with intrusions of geodes. Types F2 and F3 present a much more homogeneous surface, with colors being blue/grey and reddish/ocher. The three types correspond to the South Madrid Miocene Flint (SMM; Bustillo and Pérez-Jiménez (2005); Bustillo et al. (2012)) although their locations varied widely. Quartzite type Q1 was collected at the Jarama river (Madrid, Spain) and it is characterized by a semi-fine grained surface of light blue/grey color. Quartzites types Q2 and Q3 were also collected in fluvial deposits at Arlanzón (Burgos, Spain). Type Q2 is characterized by a grey matrix with reddish oxide bands. Type Q3 is characterized by a dark blue color and being much more fine grained than the two previous types, although some internal fissures show a minor presence of oxides. Three flakes of F1, F2, Q1 and Q2 were selected, while four flakes were selected for F3 and Q3. This makes for a total of ten flakes of flint and ten flakes of quartzite.
Te present work takes a series of steps prior to realizing any analysis and extracting data form images. These steps are: the cleaning protocol to remove contaminants, image capture, and image enhancement.
Multiple works emphasize the need of cleaning protocols to remove modern contaminants prior to analysis (Asryan and Ollé, 2020; Fernández-Marchena and Ollé, 2016; Ollé and Vergès, 2014; Pedergnana et al., 2016). Figure 1 shows the effect of manual manipulation on the texture of a surface. Because of this, a multi-step procedure based on previous studies was adopted (Pedergnana et al., 2016). This multi-step procedure included a sonic bath in 2% neutral soap (Derquim) solution during 10 to 15 minutes. This was followed by a second sonic bath in pure acetone during another 10 to 15 minutes. After each step the lithic artefacts were introduced in a water bath and finally dried using pressure air. During the cleaning protocol and microscopic analysis all artefacts were manipulated using surgical gloves.
Figure 1: Effects of manual manipulation without protection on flint surface. Left: surface photographed after the aplication of cleaning protocol and manipulation of the artefact using protection. Right: same surface after manual manipulation without protection.
Dino-Lite Edge 3.0 AM73915MZT USB microscope
-capsule
vertical stand
-pelota esa
A common problem of images obtained from USB microscopes is the lack of detail due to saturation in one of the grey level values and effects of different lightning. This saturation is often observed as a general glaze in one of the grey level values and results in a low quality image with poor detail. To increase detail and quality prior to the analysis, all images were subjected to a two-step process. First, the Fiji pluging Substract background was employed to avoid effects of different lightning. Second, the function enhance contrast was employed to desaturize the images by normalizing their histograms (Figure 2).
Figure 2: Two examples of image enhancemnt. Left: original images taken by a digital microscope. Center: images after retrieving the background to avoid effects of different lightning. Right: after normalizing the histogram to avoid saturation and increase detail.Top row: geological neocortex. Bottom row: fresh surface of experimentally knapped flake.
The present work uses three sets of statistical metrics to analyze obtained images. The first set of statistical metrics correspond to descriptive statistics (mean, standard deviation, mode, median and kurtosis). The second set of statistical measures corresponds to measures of roughness:
The implementation of these measures is done in ImageJ software (Schneider et al., 2012) through the Roughness calculation plugin (Chinga and Dougherty, 2002).
A problem exists when images are analyzed through the intensity distribution of pixels and using descriptive statistics and roughness metrics. Figure 3 illustrates the problem. Figure 3 shows a series of pattern images with similar values of intensity, but with different spatial distribution. Regarding descriptive statistics, all images present similar values of mean (127), standard deviation (125), sk (0.01) and kurtosis (-1.99). Regarding measures of roughness, all images present similar values of Rq (180), Ra (127), Rsk (1.42), and Rku (2).
Figure 3: Pattern figures presenting similar intensity distribution values, but a different spatial distribution
Figure 4: Basic representation of how the GLCM takes into consideration the spatial relationships of pixels and intensity value. The spatial relationships are considered given a reference pixel (I,J), a direction (north, east, south or west), and a distance (one pixel, two pixels)
Analysis of intensity values through the Gray Level Co-occurrence Matrix (GLCM; Haralick et al. (1973)) takes into consideration the spatial distribution of intensity values. The GLCM works in two steps (Haralick et al., 1973). First, using a given distance and direction (4), it builds a matrix which captures the relationship of intensity between pairs of pixels (reference and neighbor). Second, for every x and y it considers the co-occurrence of values, forming a new matrix. From this new matrix, a series of statistical descriptors are derived (Haralick et al., 1973).
Figure 5 shows results from the application of Haralick et al. (1973) textural features to the same pattern images.
library(tidyverse)
# GLCM of sequential data
GLCM.txt <- read.delim("Data/Pattern_GLCM.txt", sep = ",")
GLCM.txt <- GLCM.txt %>% filter(
!grepl("#", X.)) %>%
mutate(Angular.Second.Moment = as.numeric(Angular.Second.Moment),
Contrast = as.numeric(Contrast),
Correlation = as.numeric(Correlation),
Inverse.Difference.Moment = as.numeric(Inverse.Difference.Moment),
Entropy = as.numeric(Entropy)
) %>%
rename(ID = X.) %>%
select(-c(X)) %>%
group_by(ID) %>%
summarise(
ASM = mean(Angular.Second.Moment),
CONT = mean(Contrast),
CORR = mean(Correlation),
IDM = mean(Inverse.Difference.Moment),
ENT = mean(Entropy)
) %>%
mutate(
ID = substring(ID, 3)) %>%
mutate(
ID = substring(ID, 1, nchar(ID)-5))
GLCM.txt %>% pivot_longer(
cols = ASM:ENT,
names_to = "Variable",
values_to = "values") %>%
ggplot(aes("Variable", values, label = ID)) +
facet_wrap(~ Variable, scales = "free",
ncol = 2) +
xlab(NULL) +
ylab(NULL) +
geom_jitter(aes(color = ID),
width = 0.1, shape = 18, size = 3) +
ggrepel::geom_text_repel() +
ggsci::scale_color_jco() +
theme_light() +
theme(
axis.text.x = element_blank(),
strip.text = element_text(color = "black", face = "bold", size = 11),
strip.background = element_rect(fill = "white", colour = "black", size = 1),
axis.text = element_text(color = "black"),
legend.position = "none")
Figure 5: Values of pattern images using Haralick et al 1973 GLCM and second order statistics
The present work implements Haralick et al. (1973) analysis and metrics in ImageJ software (Schneider et al., 2012). Calculation of GLCM and second order statistics is done through Cabrera (2006) plugin.
Prior to the application of Haralick et al. (1973) texture metrics a proof of concept from a pilot study is undertaken. This proof of concept is aimed to two objectives. First, determine if the GLCM provides suitable variable values to discriminate between degrees of sedimentary abrasion. Second, to determine to the best interval of pixels under which the GLCM works
The proof of concept sample is composed of 23 images corresponding to geological neocortex, 56 images corresponding to flakes undergone 10 hours of sedimentary abrasion in a previous experimental study (Bustos-Pérez et al., 2019) and 42 images corresponding to fresh surfaces from experimentally knapped flakes. All flints correspond to South Madrid Miocen Flint (Bustillo and Pérez-Jiménez (2005); Bustillo et al. (2012)).
The GLCM requires to test for several distances between the reference pixel and the neighboring pixel on which the spatial relationships are established. The present work tests for five different sequences of distances (Table 1). The first sequence composes the GLCM based on pixels located at distances 1, 2 3 and 4 from the reference pixel. Sequence 2 composes the GLCM based on pixels located at distances 2, 4, 6 and 8 from the reference pixel (Table 1. All sequences of intervals are tested in the four cardinal directions (north, east, south and west). This results in 16 (four per each cardinal direction) values of Haralick et al. (1973) textural metrics. Final values are obtained by averaging the different values of each of the variables although range is not considered (Bietti, 1996).
| Sequence | Pixel Distances |
|---|---|
| 1 | 1-2-3-4 |
| 2 | 2-4-6-8 |
| 3 | 5-10-15-20 |
| 4 | 10-20-30-40 |
| 5 | 50-100-150-200 |
load("Data/PoC Data.RData")